NRG for the bosonic single-impurity Anderson model: Dynamics 



O 

(N 

Oh! 
< 

O 

(N 



i 

5-1 

CO 



I 

C 

O 
O 



> 

o 
in 

cn 
m 

o 
o 



Hyun-Jung Lee, 1 Krzysztof Byczuk, 2 and Ralf Bulla 3 

1 Asia Pacific Center for Theoretical Physics, POSTECH, Pohang, Korea 
Institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, 00-681 Warszawa, Poland 
' Institut fur Theoretische Physik, Universitat zu Koln, Koln, Germany 
(Dated: Draft: April 21, 2010) 

The bosonic single-impurity Anderson model (B-SIAM) is studied to understand the local dynam- 
ics of an atomic quantum dot (AQD) coupled to a Bose-Einstein condensation (BEC) state, which 
can be implemented to probe the entanglement and the decoherence of a macroscopic condensate. 
Our recent approach of the numerical renormalization group (NRG) calculation for the B-SIAM 
revealed a zero-temperature phase diagram, where a Mott phase with local depletion of normal 
particles is separated from a BEC phase with enhanced density of the condensate. As an extension 
of the previous work, we present the calculations of the local dynamical quantities of the B-SIAM 
which reinforce our understanding of the physics in the Mott and the BEC phases. 

PACS numbers: PACS numbers: 



I. INTRODUCTION 

The observation of the Bose-Einstein condensation 
(BEC) in ultracold, atomic gases has greatly stimulated 
research on the properties of this fascinating quantum 
state of matter X 

A particular interest lies in the controlled manipula- 
tion of the coherence and the entanglement of the BEC 
state^— which provide the basis of applications such 
as quantum computing and quantum communications. 6 
As an example, a new scheme of performing quantum 
dense coding! and teleportation^ was proposed using the 
spatial- mode entanglement of a single massive boson cou- 
pled to a BEC reservoir^. Here one considers a system 
of two coupled tightly confined potentials, each of which 
forms one of the spatial modes, A and B. In order for the 
full dense coding protocol to work, A and B have to share 
a common reference frame with which they can exchange 
particles. A BEC consisting of an indefinite number of 
particles fulfills this role 9 and the coherent control of the 
BEC state is essential to let a signal between A and B 
be phase-locked. 

On the other hand, there have been extensive studies of 
decoherence, a process of loosing quantum superpositions 
due to entanglement between a microscopic system and 
its environment. 19 The decoherence mechanism is cru- 
cial to understand the transition between quantum and 
classical system o 11 ' 12 in a sense that quantum superpo- 
sition between distinct states of macroscopic systems is 
suppressed by the decoherence process. However the en- 
vironmental effect o 10 ' 13 make it difficult to probe of the 
decoherence of macroscopic system directly. As an al- 
ternative way, one can use the coupling of a microscopic 
system, such as an atomic quantum dot, to a mesoscopic 
or macroscopic system to probe the decoherence of the 
latter. For instance, there have been several proposals 
on the single-atom-aided probe of the decoherence of a 
BECi5r.il 

In the theoretical schemes above, a BEC state is repre- 
sented as the Bose-field operator, ^ c (x) ~ /5(x) 1 / 2 e -l< ^ x - ) 



with the density p(x) and the phase 4>(x) of the conden- 
sate, of which the only available excitations at low en- 
ergies are phonons with linear dispersion. However the 
excitations of a BEC state is phononlike only for wave- 
lengths larger than the healing length £, where the heal- 
ing length £ is defined as the distance over which the 
condensate wave function grows from zero to the bulk 
value. In general, the strong collisional interaction in 
the atomic quantum dot (AQD) can locally break a BEC 
state to bring up the excitations of normal particles in- 
side of the dot. The bosonic single-impurity Anderson 
model (B-SIAM)i£ is proposed to describe the normal 
excitations in the AQD as well as the condensate part. 

Another motivation for studing the B-SIAM comes 
from a treatment of the Bose-Hubbard model within the 
dynamical mean-field theory (DMFT)^^ The DMFT 
is an exact theory in infinite spatial dimensions^ but, as 
an approximation for finite dimensional system, it was 
successful to provide comprehensive understanding about 
strongly correlated fermion systems. Recently, a new 
framework of the bosonic DMFT (B-DMFT) was pro- 
posed by Byczuk and Vollhardt^i in order to extend the 
idea of the DMFT to correlated lattice bosons and mix- 
tures of bosons and fermions on a latticed In contrast 
to the fermionic DMFT (F-DMFT), the lattice model for 
bosons is mapped into a single-impurity problem with 
two species of bath spectra, those from the condensate 
bosons and those from the normal bosons, each of which 
should be self-consistently determined. The resulting ef- 
fective bosonic impurity model has been solved by ex- 
act diagonalization (ED) method j 23 ' 24 Results have been 
presented for various phases at finite temperatures and 
compared to other theories and the experiments. 

The structure of the effective impurity model in the 
B-DMFT is reduced to the B-SIAM in the absence of 
the bath spectrum from the condensate bosons, which 
is the case in the Mott insulating (MI) phase. On the 
basis of the current work, it will be possible to perform 
NRG calculations for the B-SIAM with a self-consistently 
determined bath and investigate transitions between the 



superfluid (SF) and Mott insulating (MI) phases from 
the side of MI phase. 

The most part in this paper is devoted to discuss 
the impurity quantum phase transitions of the B-SIAM 
in terms of the local spectral density. In addition, we 
present in detail the implementation of the bosonic NRG 
for the B-SIAM to discuss various strategies to set up the 
iteration scheme for the bosonic NRG. 

The paper is organized as follows: In section [XXI we m_ 
troduce the Hamiltonian of the B-SIAM, explaining the 
differences to the spin-boson model that has been widely 
used to study the AQD coupled to a superfluid Bosc- 
Einstcin condensate^ In section IITT1 the formulation of 
the NRG for the B-SIAM is described in detail. In sec- 
tion IIV1 we discuss the impurity quantum phase transi- 
tion of the B-SIAM and explain how the Mott and the 
BEC phases are discerned in the NRG method. In sec- 
tion [V] we turn to the calculation of the local spectral 
density to discuss the different dynamical properties in 
Mott and BEC phases. Secion [VlJ is a conclusion. We 
put some technical details in appendices. 



II. MODEL HAMILTONIAN 



Va(x) 




FIG. 1: Schematic setup of an atomic quantum dot coupled to 
a superfluid atomic reservoir. The non-interacting bose par- 
ticles, denoted by the operator b k , are confined in a shallow 
trap Fb(x) and, at zero temperature, condense at the low- 
est vibrational mode to form a BEC state ^s(x) (depicted 
as a broad wave packet). The creational operator add an 
atom, (depicted by balls) in the tightly confining potential 
Va(x), where macroscopic condensation is prevented due to 
the on-site repulsion U. The atoms in Vb(x) and Vk(x) are 
coupled via a Raman transition with effective Rabi frequency 
Q. The confining potential Va(x) and Vb(x) are in all three 
directions with spherical symmetry. 



The spin-boson model has been widely used for in- 
vestigating the physical properties of an atomic quan- 
tum dot (AQD) coupled to a bosonic reser voir 1 4 i 5 1 1 5 ~ — 
In Sec. Ill A| we summarize the work by Recati et al^, 
where the particle-exchange between the AQD and the 
BEC reservoir has been discussed in terms of the spin- 
boson model. The B-SIAM model is proposed to relax 
the theoretical constrains in the spin-boson model and 
describe the density fluctuation of the coherent state orig- 
inated from the collisional interaction in the AQD. In 
Sec. IIIBl we discuss the basic set-up of the B-SIAM and 
make a comparision with the spin-boson model. 



Atomic quantum dot coupled to a superfluid 
Bose-Einstein condensate 



The particle-exchange mechanism between an AQD 
and a BEC reservoir was initially proposed by Recati 
et al. 17 The Hamiltonian is written as 



H = H f 



Hi 



Hal 



(1) 



where Hb and Ha correpond the energy of the BEC 
reservoir and the AQD, respectively. The third term 
Hab describes the Raman coupling between the AQD 
and the BEC. 

The first term Hb describes the dynamics of the BEC 
reservoir. Here the reservoir atoms are assumed to form 
a coherent matter wave, held in a shallow trapping po- 
tential Vb(x) as illustrated in Fig. [TJ The wave-function 
of the coherent state is represented as the Bose-field op- 
erator, #b(x) — p(x) 1/2 e _i ^ (x) , with the density p(x) 
and the phase 4>{x) of the condensate. 



At very low temperature, the coherent matter wave is 
regarded as superfluid Bose liquid with an equilibrium 
liquid density ps, of which the only available excitations 
are then phonons of low energy w q = « a |q| with sound 
velocity v s . In this case, the dynamics of the coherent 
matter wave is described by a hydrodynamic Hamilto- 



nian 



25 



Hb = \! dx^ PB |V0(x)| 2 



PB 



L n 2 (x) 



(2) 



where pb is the density of the superfluid fraction and 
n(x) is the density fluctuation operator II(x) = /5b (x) — 
Pb, a canonical conjugate of of the superfluid phase </>(x). 
The quadratic Hamiltonian in Eq. @ can be written in 
terms of standard phonon operators 6 q as 



H B = hv s Y,H^K 



via the following transformation, 17 



(3) 



(x) 



mv s 



ZhqVpB 



| 1/2 e lqx (& q 



n(x) = El^^l 1/2eiq - x ( fo q + fet - q )- 

2v s Vm 1 



(4) 



Here V is the sample volume. 

The second term Ha corresponds to the on-site energy 
of the AQD. The AQD is formed by trapping atoms in 
an additional tightly confining potential Va (x) as shown 
in Fig. [TJ Here one only considers the lowest vibra- 
tional mode in the AQD assuming that other higher 
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vibrational modes are off resonant due to large detun- 
ing. The collisional interaction of the atoms trapped in 
the tightly confining potential Va(x) is described by a 
coupling parameter gAA = ^TraAA^ 2 /m with scattering 
lengths qaa and atomic mass m. The strength of the 
collisional interaction between the internal states in the 
AQD and the coherent state in the BEC reservoir is given 
as gAB = iirciABfi 2 /tti. One assumes that atoms in the 
reservoir are non-interacting. Thus the on-site interac- 
tion at the AQD-site is given as, 



Ha = 



-h5 + g A B / dx|?/'h(x)| 2 / 5 B (x) 



" ( 5 ) 

where S is the detuning parameter and V'a(x) is the wave 
function of the lowest vibrational mode of the AQD. The 
on-site repulsion in the AQD is given by the parameter 
Uaa ~ 9aa/1% with I a the size of the ground state wave 
function ^(x). 

The last term Hab in Eq. (JTJ is the laser induced hy- 
bridization between particles in the AQD and the BEC 
reservoir with effective Rabi frequency fi: 



H AB = Ml I dx(* s (x)^(x) + h.c. 



(6) 



The operator V'a(x) creates an atom in the AQD and 
the operator 4 r s(x) is the annihilation operator for a 
reservoir atom at the position x. 

The Hamiltonian in Eq. (JTJ can be reduced to the 
spin-boson Hamiltonian— under the following conditions. 
First, one considers the collisional blockade limit of large 
on-site interaction Uaa, where only states with occupa- 
tion tia — and 1 in the AQD participate in the dynam- 
ics. In this case the internal state of the AQD is described 
by a pseudospin-1/2, with the spin- up or spin-down state 
corresponding to occupation by a single or by no atom 
in the AQD. Using the Pauli matrix notation, the AQD 
occupation operator b'b is then replaced by (1 + cr z )/2 
while w — > cr_|_. 

For the BEC state, one assumes that the number of 
condensate atoms inside the confinement (or the AQD) 
is much larger than 1, ns = PbI\ i- e -' ^ ne s ^ ze °f ^ ne 
spatial confinement I a is much larger than the average in- 
terparticle spacing in the BEC reservoir. Taking the long 
wave length appriximation, \<\\Ia <C 1, the phonon field 
operators in Ha and Hab are replaced by their values 
at x = 0. Further, neglecting the density fluctuations in 
the Raman coupling in Eq. ^ , the Hamiltonian Ha and 
Hab can be simplified to 



H A + H AB 



1/2 

Here A ~ fln B ' is an effective Rabi frequency. Eventu- 
ally, after a unitary transformation H = S~ 1 (Ha + Hb + 



hS 

Y 



^fn(o) J r, 



(7) 



Hab)S with S = exp{— cr z j0(O)}, the particle-exchange 
mechanism between a confined boson in AQD and a bo- 
son in the BEC reservoir can be described by the spin- 
boson Hamiltonian, 



H 



HA 



2 a °> + 

q 

-<5 + ]TA q (& q + &t) 



Here the collisional interactions and those arising from 
the coupling of the Rabi term to the condensate phase 
add coherently in the amplitudes of the phonon coupling 



Aq = 



i m ^q^ il/2 ( 9ABPB 

1 2V PB l V m v 2 s 



1 



(9) 



B. The bosonic single-impurity Anderson model 

The Hamiltonian of the B-SIAM 18 is written as 
H = £ tfb+^tfb{tfb-l) + Y,£kblb k 

k 

+ nj> t6 * + & * l 0' (io) 

k 

where b and b' are annihilation and creation operators 
obeying bosonic canonical commutation relations and 
correspond to bosons within a tight trapped potential 
Va(x), i.e. an AQD. The operators bk and b\ are an- 
nihilation and creation operators corresponding to non- 
interacting bosons confined in a shallow potential Vb( x ). 
Fig. Q] shows the schematic setup. 

The energy of the AQD is given by e and U is the local 
repulsion energy when two or more bosons occupy the dot 
system. The two parameters depend on the strength of 
the collisional interaction g a p — ATra a ph 2 /m with scat- 
tering length dap (a, j3 = A or B) and the Raman detun- 
ing 8 as discussed in Sec. Ill Al 

The third term in Eq. (|10[) is the kinetic energy of 
non-interacting bosons confined in the shallow potential 
Vb(x). Here we emphasize that the origin of the bosonic 
excitations in the B-SIAM is no more restricted to the 
phonons of the condensate wave function in the lowest 
vibrational mode in Vb(x). Instead, it involves the ex- 
cited particles to arbitrary higher vibrational modes in 
the shallow trapping potential Vg(x). The number of 
the vibrational modes in Vb(x) becomes infinite as the 
curvature of the trapping potential approaches to zero. 
In this case, the shallow trapping potential Vb(x) con- 
taining free bosons is regarded as an infinite size of a 
bosonic bath, of which the lowest vibrational mode has 
zero-energy. 

The last term in Eq. (|10l) is the laser induced hy- 
bridization between particles in the AQD and the bosonic 
bath with effective Rabi frequency f2. In analogy to the 
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fermionic SIAM the dispersion relation is determined by 
a hybridization function whose imaginary part, so called 
bath spectral function, is given by 



J(uj) = vrrj 2 ^(5(w-e fe ). 



(11) 



In the following we are interested in systems with gap- 
less bath spectral functions and in low-energy properties. 
Therefore, we use a model spectral function in the form 

J(w) = tt11 2 (1 + s) w^^VGK-w), (12) 

where Q(x) is a step- like theta function with a cut- 
off parameter w c , which yields the total spectral weight 
Jq C J(uj)duj — nil 2 . Note that the choice ui c = 1 sets 
the energy units hereafter. The exponent s characterizes 
how the bath spectral functions behave in the low-energy 
regime. 

Contrary to the spin-boson model in Ref. (17), the 
B-SIAM can consider a case where the strong Raman 
coupling f2 induces large density fluctuations around the 
AQD-site. The Raman coupling term in Eq. (fTU| im- 
poses the spatial displacement to the harmonic oscilla- 
tors in the bath, which, in consequence, increases the 
occupation of each vibrational mode. The density of the 
condensate in the lowest virbational mode increases ac- 
cordingly. Further, with Rabi coupling f2 ~ U, we go 
beyond the collisional blockade limit so that an arbitrary 
number of bosons can occupy the AQD-site to make wide 
temporal and spatial fluctuation. 

It is known that in the strong coupling regime the local 
spectrum can contain a bound or/and antibound one- 
particle states in addition to the continuum^ In this 
paper we select the coupling strength ft such that these 
extra states do not occur, which is the only restriction 
for the coupling-strength il. 

As a final remark we note that the B-SIAM Hamilto- 
nian conserves the total number of bosons. This is in 
contrast to the spin-boson model ) 27 ' 28 where the bath 
contains excited phonons, the number of which is not 
conserved. 



III. THE BOSONIC NRG 
A. Mapping onto semi-infinite chain 

In this section we describe the numerical renormaliza- 
tion group (NRG) method for conserved bosons, which 
is used to solve the B-SIAM Eq. (flOl) introduced in the 
previous section. Details of NRG for bosons are pre- 
sented in the Appendices. This method is an adoption 
of the NRG from Ref. (28) to deal with bosons with a 
conserved number of particles. 

As in the other NRG approaches^ the frequency range 
[0,w c ] of the bosonic bath spectral function J(w) is di- 
vided into intervals [aj c A~(™ +1 ), ui c A~ n ], where n = 



0,1,2,..., and A > 1 is an NRG discretization pa- 
rameter. The limit A — > 1 corresponds to the exact 
case. Within each of these intervals the spectral func- 
tion J(w c A~(" +1 ) < lj < w c A~ n ) is approximated by its 
mean value 



(uj c A- n - w c A-( ?t + 1 ))' 



(13) 



Next, following the same steps as in the spin-boson model 
in Refs. (|28ll36l) . we obtain a discretized version of the 
Hamiltonian (|10[) with new V n and e„, which are de- 
fined on a discrete frequency grid, and with new bath 
bosonic operators labeled by discrete quantum numbers 
n. Now, this discretized model is mapped onto a semi- 
infinite chai n 27 ' 28 ' 35 and we obtain the following Hamil- 
tonian 

H = etfb + —tfb(tfb-l) + V(b% + blb) 

oo oo 



The bath degrees of freedom are represented by a tight- 
binding Hamiltonian with new creation and annihilation 
operators b^} and b m , and the on-site energies e m , and 
hopping matrix elements t m between nearest neighbour 
sites. Both of them fall off exponentially, i.e. t m , e m cx 
A~ TO i 27 ' 28 ' 35 Only the first site of the semi-infinite chain, 
which is denoted by the index m = 0, is coupled to the 
impurity by the hybridization V . 

The Hamiltonian (|14p cannot be diagonalized numer- 
ically for the semi-infinite chain. Therefore, we need to 
truncate it at m — M — 2, which corresponds to tak- 
ing M sites, including the impurity-site, in the chain. 
Since the Hamiltonian parameters t m , e m decay expo- 
nentially with m, this truncation is justified at large M. 
The Hamiltonian diagonalized numerically has the form 

H M = etfb+^-tfb{tfb-l) + V(b% + P b) (15) 



M-2 



Af-3 



+ emblem + Y tmibtnbm+l + b\ n+1 b m ) . 



m=0 



m— 



The Hamiltonian (|15|) commutes with the number oper- 
ator 



N 



M 



¥b- 



M-2 

E 

m— 



bL b r . 



(16) 



Hence, the eigenstates of Hm are also the eigenstates of 
Nm, so they are labeled by the corresponding quantum 
number N. The Hilbert space of all states with the same 
N is denoted by T-Ln- The dimension of each Hilbert 
space T-Lm with a given M is 



V 



N 



(M-l + N)\ 
(M — 1)! N\ 



(17) 
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Unfortunately, for large N and M the Hilbert space di- 
mension is so large that direct diagonalization methods 
are not efficient. Therefore, the Hamiltonian (fT5|) is di- 
agonalized iteratively as is discussed next. 



B. Iterative Diagonalization 

At the beginning for small M and N such that the 
Hilbert space dimension Djsr is less than typically few 
thousands, which depends on the computing facility, we 
perform exact diagonalization of the Hamiltonian (fT5)) 
for a given M and all possible N such that 



N = 0,1,2,. ..,N n 



(18) 



where A max is a cutoff for a number of particles. The 
truncation of the possible particle numbers, which is an 
approximation, is a necessary to make a computation fea- 
sible. As we will see later if the cutoff A max is large 
enough then it does not affect obtained results. 

Having diagonalized the Hamiltonian Hm for a given 
M we increase the system size by adding one more site to 
the chain. Then we diagonalize the Hamiltonian Hm+i 
which has the form 



Hm+i — Hm + £iu-ib\ I _ 1 b 



M-l 



(19) 



tM-2(b\. I _ 2 bM-l + b M _ x bM-2) 



If it turns out that the dimension of the Hilbert space is 
too large now, we need to construct an effective represen- 
tation of the low energy eigenstates while M increases. 
This is done iteratively as is described below. 

We keep the dimension of the Hilbert space constant 
by taking only the low energy eigenstates. However, to 
be able to make a direct comparison of the spectra while 
M increases we need to scale the M + 1 site Hamiltonian 
as follows 



H 



M+l 



AH 



M 



A 



Af-l 



£M-ib M _ 1 bM-i- 



(20) 



tM-2(b M _ 2 bM-l + b M _ x bM-2) , 



where we keep the same symbol for the Hamiltonian. All 
eigenvalues of Hm for all < N < N max are sorted in 
an ascending way, and the N s eigenstates | N, t^)m with 
the lowest eigenvalues are used in diagonalizing Hm+i- 
Explicitly, we take into account such states that 

H M \N,r N ) M = E rN M(N)\N,r N ) M , (21) 



with rjv = 1, nY" ', where ni N ^ is the number of N- 
particle states with the lowest eigenvalues E Tn> m(N) in 
each Hilbert space T-Ln^ The dimension of the Hilbert 
space N s is given by the summation of n 



(JV) 



N s = r> 

N=0 



(N) 



(22) 



and optimized to perform the computation feasible^ 

In the Hilbert space of the Hm+i Hamiltonian, the 
A-particle states are given by 

{\N, R) M +i} = {\N - k, r N - k )M ® |fc)}fe=o,...,JV, (23) 



where 



\k) = v M zl J |0), 



(24) 



is the fc-particle state on the M — 1 site in the chain, 
and |0) is an empty (vacuum) state on this last site. In 
Eq. (|23[) the quantum number R = r^-k means the quan- 
tum number of the Hm Hamiltonian with N—k particles. 
The numbers R are not the quantum numbers labeling 
the eigenstates of the Hamiltonian Hm+i- This is due 
to the fact that the new Hamiltonian Hm+i does not 
commute with the total number of particles Nm of the 
previous system with the Hamiltonian Hm, i-e. we can 
check that [Hm+i, Nm] ^ 0, where Nm is defined in (flUj) . 
In order to find eigenstates of Hm+i in a basis (f2"3"|) we 
construct the Hamiltonian matrix elements 



H{R;R') = M +i (N,R\H M +i\N,R') 



M+l, 



(25) 



and diagonalize this matrix obtaining a set of eigenvalues 
and eigenstates 



\N,u n )m+i = ^U N (u N ;R)\N,R) 



M+l 



(26) 



where Un (wjv ; R) is an orthogonal matrix, and wv are 
new quantum numbers labeling an N particle eigenstate 
of Hm+i with eigenvalue E UNt M+i{N). The procedure 
described from Eqs. (|2"3")l to (f2"Bf is repeated for all N = 
0, 1, 2, N max . 

In the next iteration step we extend the system by 
adding one more site to the chain and use the eigenstates 
(|2"6")l of Hm+i to construct a basis of the new Hamilto- 
nian in a way analogous to (|23[) . Repeating the same 
procedure as described between Eqs. (|2T1) to (|2"6")l we ob- 
tain new eigenstates and eigenvalues of the larger sys- 
tem. Further details on the iterative diagonalization is 
presented in Appendix [A"l 

We proceed iterative diagonalizations until the many 
particle spectra approach the trivial fixed point of the 
non-interacting bosonic bath. The low-energy spectrum 
of Mott and BEC phases and the structure of the fixed 
points are presented in Sec. IIVB1 



IV. ZERO-TEMPERATURE PHASE-DIAGRAM 

A. Overview 

The zero-temperature phase diagram in Fig. [5] is cal- 
culated for fixed U = 0.5uj c with the parameter space 
spanned by the dimensionless coupling constant a = 



6 
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FIG. 2: Zero-temperature phase diagram of the B-SIAM for 
bath exponent s = 0.4 and fixed impurity Coulomb inter- 
action U = 0.5uu c - The different symbols denote the phase 
boundaries between Mott phases and the BEC phase. The 
Mott phases are labeled by the number of the impurity- 
quasiparticle, ni mp . Only the Mott phases with n; mp < 4 
are shown. The NRG parameters are A = 2.0, Ni, = 10, and 
N a = 100. 



(a) Reference frame without impurity 
Vb(x) 




FIG. 3: (a) Reference frame: A non-interacting BEC state, 
$s(x), is confined in a shallow trapping potential Vb(x). 
(b) Mott phase: The impurity-quasiparticle consists of an in- 
teger number of depleted particles, (depicted as balls) , which 
are tightly trapped in Va(x). The other bosons contained in 
Vb (x) still form a BEC cloud but the local density of the con- 
densate vanishes in the vicinity of the AQD. (c) BEC phase: 
The impurity-quasiparticle forms a part of a BEC state to 
enhance the density of the condensate around the AQD. The 
confining potential Va (x) and Vb (x) are in all three directions 
with spherical symmetry. 

■^ii^-Sl 2 and the impurity energy e. We choose s = 0.4 
as the exponent of the power law in J(oj) in Eq. (fT2|) . A 
similar phase diagram for different bath exponent s = 0.6 
has been presented in Ref. (|18f ). The phase diagram is 
characterized by a sequence of lobes. We use the termi- 
nology 11 Mott phase" for the inside of the lobes and "BEC 
phase" for the region outside of the lobes. 

The Mott and the BEC phases are distinguished by 



a hybridized state that is formed around the AQD as 
illustrated in Fig. [3] Fig. |3J-(a) shows a BEC state of 
a non-interacting bosonic bath, where all existing parti- 
cles occupy the lowest vibrational mode of the shallow 
potential Vb(x). In the presence of the AQD, however, 
particles around the AQD can be either completely de- 
pleted (Fig.[3]-(b)) or even more concentrated toward the 
local site (Fig. 0]-(c)). We call the collective excitation 
around the AQD as impurity-quasiparticle. 

In the Mott phase, the impurity-quasiparticle consist 
of an integer number of depleted particles, (depicted as 
balls in Fig. G2-(b)), which are tightly trapped in Va(x). 
The number of the depleted particles is used to label the 
different Mott phases in the phase diagram in Fig. [5] The 
other bosons contained in Vb(x) still form a BEC cloud 
but the local density of the condensate vanishes in the 
vicinity of the AQD. 

In the BEC phase (Fig. |3J-(c)), the impurity- 
quasiparticle forms a part of a BEC state to enhance the 
density of the condensate around the AQD. The enhance- 
ment of the condensate-density is due to the strong Ra- 
man coupling and the deep attractive potential e < 
of the AQD. 

Numerical evidences for our assertions are presented 
in the rest part of the paper. In Sec. IIVB1 we look 
into the contribution of the impurity-quasiparticle to the 
ground state energy. In Sec. [V] the calculation of the 
local Greens function is presented to show the local dy- 
namics of normal and condensate particles. 



B. Impurity contribution to the ground state 
energy 

The ground state energy of a non-interacting bosonic 
bath is zero since all existing particles occupy the low- 
est vibrational mode with zero-energy. An impurity site 
with repulsive interaction U, however, can deplete some 
particles from the zero-energy mode and shift the ground 
state energy to be finite. In general, the non-zero ground 
state energy depends on the total number of particles 
(N) in the system. 

Fig.|4]and Fig.[5]show the A^-dependence of the ground 
state energy E$^m{N) in Mott and BEC phases, respec- 
tively. The different curves are the results from differ- 
ent size (M) of the systems. The ground state energy 
Eq,m{N) decreases until the configuration around the 
AQD, (i.e. impurity-quasiparticle) is optimized. The oc- 
cupation at the minimum point is denoted by N*. 

The minimum ground state energy Eq^m(N) at N — 
N* is plotted as a function of the system size M in Fig. [6] 
The minimum ground state energy Eq^m{N) at N = N* 
converges in the limit M — > oo. Once the system con- 
verges into the large M limit, all ground states for differ- 
ent N become degenerate. Indeed, Fig.|4]and Fig.[5]show 
that the ground state energy Eq^m{N) becomes almost 
independent of N already for M = 9. 

In the thermodynamic limit (N —¥ oo, M —¥ oo), the 
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FIG. 4: Ground state energy Eo,m(N) vs. JV calculated for 
s = 0.6, U = 0.1, e = -0.36 and V = 0.01. The parameters 
are chosen inside of the Mott lobe labeled by 4 (Mott phase 4) . 
The inset shows the position of the minimum point N = N* 
as a function of M. The NRG parameters used are A = 1.5 
and N. = 1000. 
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FIG. 5: Ground state energy Eo,m(N) vs. TV calculated for 
s = 0.6, U = 0.1, e = -0.05 and' V = 0.4. The parameters 
are chosen outside of the Mott lobes (BEC phase). The in- 
set shows the position of the minimum point N — N* as a 
function of M. The NRG parameters used are A = 1.5 and 
N s = 1000. 




FIG. 7: The lowest lying many-particle levels E n ,M^- ver- 
sus iteration number M for parameters s — 0.7, V = 0.01, 
U = 0.5, and e = -1.2 (Mott phase). The NRG parameters 
used are A = 1.25, N s = 3000 and N max = 15. 




— V=0.40, 11=0.1, e=-0.05 (BEC phase) 
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FIG. 6: The minimum of the ground state energy Eo,m{N) 
at N = N* as a function of M. The solid and dashed lines 
correspond to the BEC and the Mott phases, respectively. 
The bath exponent is fixed to s = 0.6. The NRG parameters 
used are A = 1.25, N max = 40 (8) and N s = 8000 (600) for a 
BEC (Mott) phase. 



FIG. 8: Flow diagram of the lowest lying many-particle levels 
£n,MA M_1 versus iteration number M for parameters s — 
0.7, V = 0.4, U = 0.1 and e = -0.05. The NRG parameters 
used are A = 1.25, N s = 8000 and N max = 40. 

result of adding (or removing) one particle is to convert 
a state of a system of N particles into the same state of 
a system of N ± 1 particles: 

lim lim |iV±l,0) M = lim lim |iV,0V. (27) 

N^oo M— s-oo N—>oo M— >oo 

Here \N, 0)m is the TV-particle ground states of Hm ■ This 
is the case of a condensate consisting of a macroscopic 
number of particles, i.e. a coherent statei--i 

The degenerate feature of the ground states in Eq. (|2"7|) 
is extended to the low lying excited states when the 
many-particle spectrum reaches a fixed point. 

Fig. [7] shows the energy flow of the lowest lying many- 
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FIG. 9: The quasiparticle-occupation N* versus the system- 
size M. The data obtained for parameter parameters s — 0.7, 
V = 0.4, U = 0.1, and e = -0.05 (BEC phase). The NRG 
parameters are A = 1.25, N max = 40 and N s — 8000. 



particle levels E n ^M{N) versus iteration number M . The 
parameters V, U and e are chosen for the system to flow 
into a Mott phase. Three pannels show the iV-particle 
eigenstates for TV = 9, 10, and 11. The eigenstates in the 
three figures flow into the same fixed point, which is a 
trivial fixed point of a non-interacting bosonic bath. It 
means that the dynamics of the AQD, i.e. the impurity- 
quasiparticle, is suppressed in this energy-scale so that 
the low-lying excitations show the dynamics of the non- 
interacting bosons that locate far from the AQD-site. 

Fig. [5] shows the lowest lying many-particle levels in a 
BEC phase. Three pannels show the iV-particle eigen- 
states for TV = 26, 29, and N = 32, which flow into the 
same strong-coupling fixed point. The level-spacing in 
the strong-coupling fixed point is different from the one 
in the non-interacting fixed point - the reason is not clear 
yet. 

As a last remark, we mention the conditions for nu- 
merical convergence. We see that the energy-levels start 
to deviate from the strong-coupling fixed point around 
at the iterative step M = 20. The upturn (deviation 
from the fixed point) appears if the number of parti- 
cles iV is not large enough compared to N*(M). The 
N* increases with increasing M (see Fig. ^ and reaches 
the value N* ~ 30 at the iteration M = 20. The N- 
particle eigenstates flows into the same strong-coupling 
fixed point only if TV is larger than N* . 

The quick and the slow convergence in Mott and BEC 
phases respectively can be interpreted as following. The 
system-size M corresponds to the number of vibrational 
modes that are taken into account in Hm, i.e. the larger 
system involves more vibrational modes with small en- 
ergy. From this we can conclude that the impurity- 
quasiparticle in a Mott phase consists of the depleted par- 
ticles occupying the higher vibrational modes in Vb(x), 
which can be described by a relatively small system size. 
In a BEC phase, however, the impurity-quasiparticle is 
a part of a condensate which consists of a macroscopic 
number of particles with almost zero-energy. Thus one 
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FIG. 10: (a) The local spectral density of the B-SIAM for 
bath exponent s = 0.6 and fixed impurity Coulomb interac- 
tion U = 0.5w c , the onsite impurity energy e = —0.7, and the 
hybridization V = 0.15 (Mott phase 2). The NRG parameters 
are A = 1.25, N max = 3, and N s = 1000. (b) The (positive) 
low-frequency part of A(uj) (solid line) in log-log scale. The 
dashed line is a guide line for eyes showing a power-law be- 
havior (oc oj s , s — 0.6). (c) The low- frequency part of A(oj) 
in linear scale. A(u) vanishes at ui — 0. 



needs a large value of N and M to properly describe the 
condensate. 



V. LOCAL DYNAMICS AT ZERO 
TEMPERATURE 

The local Green's function of the impurity model is 
defined as 



G{z) 



1 



dte izt {[b{t),tf]}. 



(28) 



where b^' is an annihilation (creation) operator for the 
impurity. The local spectral density A(uj) is the imagi- 
nary part of the local Green's function, 



A{uj) = --%G(uj + i5). 



(29) 



The local spectral density in a Mott phase (Fig. [TU)) 
shows two quasiparticle peaks that are separated by a 
gap, Ag a p ~ 0.2. A sharp peak at ui ~ —0.2 is a signal of 
hole-excitation in the AQD. The particles trapped in the 
AQD show no resonance with the reservoir as if they are 
isolated from it. In fact, most of particles in the reservoir 
are immobile since they are condensed at zero-energy and 
make no resonance with the particles in the AQD. 

The local occupation at the AQD-site can be obtained 
by integrating the spectral weight below the chemical po- 
tential fi = 0, 



ni oc (T = 0) 



f B E{u)A(u))dio 



T=0 



1.8875 



(30) 
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FIG. 11: (a) The local spectral density of the B-SIAM for bath 
exponent s — 0.6 and fixed impurity Coulomb interaction 
U — O.lcjc, the onsite impurity energy e = —0.05, and V = 0.3 
(BEC phase). The NRG parameters are A = 1.25, N max = 
40, and N„ = 5000. (b) The low- frequency part of |A(u)| as a 
function of \cj\ in log-log scale, where upper and lower curves 
correspond to the positive and negative spectral density in the 
first and third quadrant, respectively. The dashed line in the 
inset is a guide line for eyes showing a power-law behavior 
(oc lj~ s , s — 0.6). (c) The low-frequency part of A(co) in 
linear scale. A(ui) shows two 5-peaks at cj ~ ±0.0001. The 
position of two peaks approach u — in thermodynamic limit 
M — > oo as seen in Fig. 1131 



where the Bose-Einstein distribution function is 
given as a step function at zero temperature, 



lim /bb(w) = lim 



1 



1 



= -e(-w) 



(31) 



with 6(w) = 1 for u) > and Q(uj) = for u> < . 

Creating a particle at the AQD-site gives a broad peak 
at positive frequency. There is no feature at uj = 
(Fig. HDJ-(c)) indicating that the BEC is locally forbid- 
den around the AQD-site. The A(ui) vanishes at uj = 
with a power-law behavior, 



A(uj) kw s , uj > 0, 



(32) 



which is the same for the bath spectral function J(oj) 
in Eq. (|12p. The power-law behavior for various bath 
exponents s is shown in Fig. 1121 (a). 

Fig. [TT] shows the spectral density in the BEC phase. 
The spectral density in the BEC phase (Fig. [TT]) diverges 
at uj = 0, 



A{uj) oc sgn(w)|w| 



(33) 



where the power-law corresponds to the inverse of the 
bath spectral density J(oj), see Fig. [UJ-(b). The diver- 
gence of A(u) occurs if a hybridized state is pinned at the 
gapless point of the spectral function J(oj). To discuss 



more details, let us look into the local Green's function 
G{z) written as 

G(z) = (z-e- Z(z))- 1 (34) 

where e is the energy of the impurity level (with operator 
b(") and is the total self-energy of the impurity 

model. The imaginary part of the Green's function in 
Eq. ([34]) is given as 



Z[G(z)} = 



9[E(*)] 



S(z)])2 + (3[S(z)])2 



(35) 



(3? [z- 

with z = uj + i0 + . The actual calculation of is in 
progress and will be presented in our subsequent paper. 
Here we assume that the imaginary part of the self-energy 
[S(z)] follows a power-law behavior with the same ex- 
ponent as the bath spectral function J(ui) oc uj s : 

3[E(z)]ocw s . (36) 

The singular behavior of the local spectral density A(lu) 
in Eq. (|33[) can appear when the impurity bound state 
occurs at uj = 0: 



3? [u 



E(ui +i0 + )] = at uj = 0. 



(37) 



The imaginary part of the self-energy shows a power-law 
behavior as assumed in Eq. (|3"6"|) . In the case, the 3 [G(z)] 
becomes inverse-proportional to 3 [£(z)], 



3 [G{z)\ oc 



1 



9 [£(*)] 



(38) 



If the impurity bound state occurs below the chemical 
potential, the first term in the numerator in Eq. (|35|) is 
non-zero at uj = 0, which makes 3 [G(z)] proportional to 
3 [S(z)] around the gapless point uj — 0, 



9 [G{z)\ oc 3 [E(z)] oc uj s 



(39) 



A similar feature of A(w) is observed in the pseudo-gap 
Anderson model ) 30 i 40 where a Kondo bound state ap- 
pears at the gapless Fermi level. 

The singular behavior of A(uj) for various bath expo- 
nents s is shown in Fig. H2l-(b). 

Another interesting feature in the BEC phase is the 
finite spectral weight at uj = as shown in Fig. [TTI-(c). 
Fig. ITH-(c) shows two peaks at small frequency ujq± m 
±0.0001 with opposite sign of spectral weight. In the 
limit M — > oc, the position of both peaks approaches to 
zero (oj = 0) and the amplitude |7o±| converges to the 
same value (Fig. [13]) . The finite spectral weight at uj = 
indicates the existence of the condensate particles in the 
AQD-site. 

The local occupation at the AQD-site can be obtained 
from integrating the spectral weight below the chemical 
potential fi = 0, 



ni oc (T = 0) 



f BE (uj)A(uj)duj 



T=0 



7o- 
4.6 



h lim 

e->0- 

0.22 



A(uj)du 



(40) 
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FIG. 12: (a) The low-frequency behavior of A(uj) for U = 
0.5lu c , e — —0.7, V = 0.15 (Mott phase 2) and for various 
bath-exponent s = 0.5, 0.6, and 0.7. (b) The low-frequency 
behavior of A(lj) for U = 0.1w c , e = -0.05, V = 0.3 (BEC 
phase) and for various bath-exponent s = 0.5, 0.6, and 0.7. 
The NRG parameters are A = 1.25, N max = 40, and N s = 
5000. The dashed lines are guide lines for eyes to show the 
power-law behavior. The NRG parameters are A = 1.25, 
N max = 3, and N s = 1000. 




M 



FIG. 13: The position (\u)o± |) and the amplitude (|7o± |) of the 
two peaks in Fig. 1111 (c) depending on the size of system M. 
The indices 0+ (0—) denote the peak at positive (negative) 
frequency, respectively. 



where the Bose-Einstein distribution function at 
zero temperature is given in equation (I31[) . The first 
term (= 4.6) in equation (|40|) is the contribution of the 
condensate particles whereas the second term(= 0.22) is 
the contribution of particles that are depleted from the 
condensate. 



VI. CONCLUSION 

The bosonic single-impurity Anderson model is studied 
to understand the local dynamics of an atomic quantum 
dot (AQD) coupled to a BEC state. The major result 
presented in this paper is the calculation of the impu- 
rity Green function but, in addition, considerable space 
is devoted to refine the description of the Mott and the 
BEC phases. The local collisional interaction, dominant 
over the Raman coupling, depletes the particles around 
the AQD out of the condensate (Mott phase). Other- 
wise, the Raman transition makes the density of the BEC 
state even more concentrated toward the local site (BEC 
phase). The AQD can share a coherent phase of the 
macroscopic condensate only in the BEC phase and can 
be used to probe the decoherence of the BEC stated - — 

The scheme for the quantum dense coding protocol^ 
requires two separate AQDs, both of which are coupled 
to the same BEC state. In Ref. (4), it is assumed that 
a signal between the two AQDs is phase-locked through 
a BEC state with uniform density and phase. However 
the phase preserved in each AQD can depend on the po- 
sition of the dots when the AQDs make the BEC state 
non-uniform. In this case, the spatial fluctuation of a 
BEC cloud in the presence of two AQDs deserves of fur- 
ther research, for which a recent extension of the NRG 
technique, computing spatial correlation function for the 
Kondo screening cloud^ is also applicable. 
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Appendix A: Details about the Iterative 
Diagonalization 

Now we obtain the matrix elements in Eq. (|25[) . 

H{R;R') =m+i (N,R\H m+ i\N,R')m + i, (Al) 

where the TV-particle states \N,R)m+i is defined in 
Eq. (jUD. 

It is straightforward to demonstrate that the diagonal 
matrix elements of Hm+i arc 

H(R; R) = E RM {N - k) + k e M -i- (A2) 
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The only non- vanishing off-diagonal elements of Hm+i 
are given by 

H(R';R) 

= 4',fc-iiAf- 2 Vfc M (N - k', R'\\P M _ 2 \\N ~ k, R) M 

+Sk>.k+it M -2\ / k + 1 m(N - k', R'\\b M -2\\N - k, R) M , 

(A3) 

where (H&M-2II) are ^ ne invariant matrix elements. 

In obtaining Eq. (|A2I) , we have made use of the follow- 
ing results, 

M + i{N , ,R'\\b M -i\\N,R)M+i = S k , ik - 1 Vk (A4) 

and 

M+l (N'^R'fblj^WN.R) M+l — <5fc',fc+l VkTT. (A5) 

which follow from the definition of the basis set in 
Eq. ([23]). 

From Eq. (|A2|) and Eq. (|A3|) . it is clear that we can 
set up the matrix of H (R; R') starting with the knowl- 
edge of the previous iterative step such as the eigenenergy 
En t M{N — k) and the matrix elements 

M (N -k- l,r / ||5 Ar _a||JV - k,r) M (A6) 

for k = 0,...,N. 

The actual iteration upon entering the stage (M + 1) 
would proceed as follows. We first start with the lowest 
allowed value of Nm+i(= 0), and then increase it in steps 
of 1. Within a given /Cjv subspace, we construct the 
matrix 

H(R; R') =m+i (N, R\H M +i\N, R') M +x- (A7) 
Diagonalization of this matrix gives a set of eigenstates 

\N,u N ) M+ x =J2U N (u N ]R)\N,R) M +i (A8) 

R 

where Un will be an orthogonal matrix. The diagonaliza- 
tion means no more than the knowledge of m+i(N) 
and J7/v(wjv; R). After completing the diagonalization 
for one N, we proceed up, increasing TV in steps of 1. 
In order to go to the next iteration we need to calcu- 
late m+i(N — l,uj'\\bM-i\\N,uj)M+i- Using the results 
in Eq. (|A4[) , it is easy to verify that 

m+i(N - l,u / N _ 1 \\b M -i\\N 1 u! N )M+i 
= J2 u N-^N-i\R)UN(u} N ;R)Vk (A9) 

R 

where k is the number of particles on the M — 1 site in 
the chain as given in the Eq. (f2"3"|). 
Appendix B: Calculation of local spectral density 

The NRG method uses a discretized version of the An- 
derson model in a semi-infinite chain form in Eq. (1151) . 



The resulting spectral functions will therefore be given 
as a set of discrete (5-peaks. For example, the spectral 
representations of the one-particle Green's function G(z) 
is 

A(u) - --$SG(u) 

= E E K^r||&tj|AT',r'>| 2 exp{-^(iV 5 r)} 

N,r N',r' 

x5(u} — E(N, r) +E{N',r')) 
"E E K^V, r| |6| |A^', r') | 2 exp {-^(A^, r)} 

N,r N'.r' 

x5(lj + E(N, r) — E{N' , r')) 

(Bl) 

Here \N,r) and E(N,r) are the abbreviation of \N,t)m 
smdE M (N,r) in Eq. flip. 

As a practical matter, however, calculating the states 
of Hm for large N is hard to deal with because the num- 
ber of iV-particle states of Hm explodes in combinatorial 
way as shown in Eq. (1171) . Thus we introduce cut-off, 

N maac 

E - E ■ (B2) 

n n=o 

The value of N max has to be larger than the minimum 
point of the ground state energy at N = N* . 

At zero temperature, the ensemble average in 
Eq. (|B1[) is replaced to the ground expectation value 
(JV*,0|...|JV*,0) : 

A{u)) T=Q = --3G(w) T=0 

7T 

= ^|<7V* + l,r|6t|7V*, )| 2 

r 

x5(u - E(N* + 1, r) + E(N*,0)) 
-^|(7V*-l,r|6|A*,0)| 2 

r 

xS{lu + E(N* - l,r) - E{N*,0j). 

(B3) 

The matrix elements (N, r\\w\\N', r') and the energies 
E(N, r) are calculated in the NRG method. The resulting 
spectral function, as a set of <5-functions at frequencies w n 
with weights g n , are broadened on a logarithmic scale as 

e -b "/ 4 (Inui — lnw„) 2 . 

g n 5(u) - u n ) -> 9n7 7= exp -ny . (B4) 

b n u n ^/TT b n 

In our calculations, the width b n is chosen as b inde- 
pendent of n and the typical values we use are in the 
range 0.01 < b < 0.1. A (5-peak in Fig. flUKa) is an in- 
trinsic (5-peak without any resonance, for which we use a 
value, b n = 0.0001. 
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A(e) oc£ r e(|e/D- 1|). 



(B5) 



Here D is a band- width of the host metal and 6(e) is 
a step-like theta function. In the strong-coupling phase 
(particle- hole symmetric case), a Kondo bound states 
pinned at the Fermi level e = 0, where A(e) becomes gap- 
less. The local spectral density (u>) diverges at lj = and 
the low-frequency behavior is given as 3 -^ 



A(u) 



(B6) 



